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I.  \ Introduction 

^ The  objective  of  this  investigation  is  to  contribute  to 
numerical  techniques  of  interest  in  system  theory  in  general 
and  aerospace  engineering  in  particular.  Typical  areas  of 
matehmatical  research  ares  Urf  solution  of  nonlinear,  algebraic 
or  trascendental , equations/  £bf  solution  of  nonlinear,  ordi- 
nary differential  equations/  (c)  solution  of  linear  and  non- 
linear partial  differential  equations ) fdT*  mathematical  pro- 
gramming problems;  and  (erf  optimal  control  problems.  Typical 
areas  of  application  are:  fear)'*'' numerical  solution  of  optimum 
structural  problems,’  (-iirf  numerical  solution  of  optimum  flight 
trajectories/  and  (x±±T  numerical  solution  of  optimum  aero- 
dynamic shapes.  Algorithms  for  digital  computer  usage  are 
being  developed  and  tested  on  IBM  370/155. 
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V.  Abstracts  of  Reports  Completed 

1.  CLOUTIER,  J.R.,  MOHANTY,  B.P.,  and  MIELE,  A.,  Sequential 
Conjugate  Gradient-Restoration  Algorithm  for  Optimal 
Control  Problems  with  Nondifferential  Constraints, 

Part  1,  Theory,  Rice  University,  Aero-Astronautics 
Report  No.  126,  1977. 

Abstract.  A sequential  conjugate  gradient-restoration 
algorithm  is  developed  in  order  to  solve  optimal  control  pro- 
blems involving  a functional  subject  to  differential  con- 
straints, nondifferential  constraints,  and  terminal  constraints. 
The  algorithm  is  composed  of  a sequence  of  cycles,  each  cycle 
consisting  of  two  phases,  a conjugate  gradient  phase  and  a 
restoration  phase. 

The  conjugate  gradient  phase  involves  a single  iteration 
and  is  designed  to  decrease  the  value  of  the  functional  while 
satisfying  the  constraints  to  first  order.  During  this  iter- 
ation, the  first  variation  of  the  functional  is  minimized,  sub- 
ject to  the  linearized  constraints.  The  minimization  is  per- 
formed over  the  class  of  variations  of  the  control  and  the  pa- 
rameter which  are  equidistant  from  some  constant  multiple  of 
the  corresponding  variations  of  the  previous  conjugate  gradient 
phase.  For  the  special  case  of  a quadratic  functional  subject 
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to  linear  constraints,  various  orthogonality  and  conjugacy 
conditions  hold. 

The  restoration  phase  involves  one  or  more  iterations 
and  is  designed  to  restore  the  constraints  to  a predetermined 
accuracy,  while  the  norm  of  the  variations  of  the  control  and 
the  parameter  is  minimized,  subject  to  the  linearized  cons- 
traints. The  restoration  phase  is  terminated  whenever  the 
norm  of  the  constraint  error  is  less  than  some  predetermined 
tolerance. 

The  sequential  conjugate  gradient-restoration  algorithm 
is  characterized  by  two  main  properties.  First,  at  the  end 
of  each  conjugate  gradient-restoration  cycle,  the  trajectory 
satisfies  the  constraints  to  a given  accuracy;  thus,  a sequence 
of  feasible  suboptimal  solutions  is  produced.  Second,  the  con- 
jugate gradient  stepsize  and  the  restoration  stepsize  can  be 
chosen  so  that  the  restoration  phase  preserves  the  descent 
property  of  the  conjugate  gradient  phase;  thus,  the  value  of 
the  functional  at  the  end  of  any  cycle  is  smaller  than  the 
value  of  the  functional  at  the  beginning  of  that  cycle.  Of 
course,  restarting  the  algorithm  might  be  occasionally  neces- 
sary. 

To  facilitate  numerical  integrations,  the  interval  of 
integration  is  normalized  to  unit  length.  Variable- time 
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terminal  conditions  are  transformed  into  fixed-time  terminal 
conditions.  Then,  the  actual  time  at  which  the  terminal 
boundary  is  reached  becomes  a component  of  a vector  parameter 
being  optimized. 

Convergence  is  attained  whenever  both  the  norm  of  the 
constraint  error  and  the  norm  of  the  error  in  the  optimality 
conditions  are  less  than  some  predetermined  tolerances.  Seve- 
ral numerical  examples  illustrating  the  theory  in  this  paper 
are  given  in  Part  2. 

Key  Words.  Optimal  control,  gradient  methods,  conjugate- 
gradient  methods,  numerical  methods,  computing  methods, 
gradient-restoration  algorithms,  sequential  gradient-restoration 
algorithms,  sequential  conjugate  gradient-restoration  algorithms, 
nondifferential  constraints. 


2. 


CLOUTIER,  J.R., 


MOHANTY,  B.P.,  and  MIELE,  A.,  Sequential 


Conjugate  Gradient-Restoration  Algorithm  for  Optimal 
Control  Problems  with  Nondifferential  Constraints, 

Part  2,  Examples,  Rice  University,  Aero-Astronautics 
Report  No.  127,  1977. 

Abstract.  In  Ref.  1,  Cloutier,  Mohanty,  and  Miele  de- 
veloped the  sequential  conjugate  gradient-restoration  algorithm 
for  minimizing  a functional  subject  to  differential  constraints, 
nondifferential  constraints,  and  terminal  constraints.  In  this 
report,  several  numerical  examples  are  presented,  some  pertain- 
ing to  a quadratic  functional  subject  to  linear  constraints  and 
some  pertaining  to  a nonquadratic  functional  subject  to  non- 
linear constraints.  These  examples  demonstrate  the  feasibility 
as  well  as  the  rapid  convergence  characteristics  of  the  sequen- 
tial conjugate  gradient-restoration  algorithm. 

Key  Words.  Optimal  control,  gradient  methods,  conjugate- 
gradient  methods,  numerical  methods,  computing  methods, 
gradient-restoration  algorithms,  sequential  gradient-restoration 
algorithms,  sequential  conjugate  gradient-restoration  algori- 
thms, nondifferential  constraints. 


PR-11 


3.  MIELE,  A.,  GONZALEZ,  S.,  and  WU,  A.K. , On  Testing 

Algorithms  for  Mathematical  Programming  Problems,  Rice 

University,  Aero-Astronautics  Report  No.  134,  1976. 

Abstract.  This  paper  considers  the  comparative  evalua- 
tion of  algorithms  for  mathematical  programming  problems.  It 
is  concerned  with  the  measurement  of  computational  speed  and 
examines  critically  the  concept  of  equivalent  number  of  fun- 
ction evaluations  Ng.  Does  this  quantity  constitute  a fair 
way  of  comparing  different  algorithms? 

The  answer  to  the  above  question  depends  strongly  on 
whether  or  not  analytical  expressions  for  the  components  of 
the  gradient  and  the  elements  of  the  Hessian  matrix  are 
available.  It  also  depends  on  the  relative  importance  of  the 
computational  effort  associated  with  algorithmic  operations 
viv-a-vis  the  computational  effort  associated  with  function 
evalauations. 

Both  theoretical  considerations  and  extensive  numerical 
examples  carried  out  in  conjunction  with  the  Fletcher-Reeves 
algorithm,  the  Davidon-Fletcher-Powell  algorithm,  and  the 
quasilinearization  algorithm  suggest  the  following:  the  Ng 
concept,  while  accurate  in  some  cases,  has  drawbacks  in  other 
cases;  indeed,  it  might  lead  to  a distorted  view  of  the  rela- 
tive importance  of  an  algorithm  with  respect  to  another. 

The  above  distortion  can  be  corrected  through  the  in- 
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troduction  of  a more  general  parameter  Nq.  This  generalized 
parameter  is  constructed  so  as  to  reflect  accurately  the 
computational  effort  associated  with  function  evaluations  and 
algorithmic  operations. 

From  the  analyses  performed  and  the  results  obtained, 
it  is  inferred  that,  due  to  the  weaknesses  of  the  Ng  concept, 
the  use  of  the  Ng  concept  is  advisable.  In  effect,  this  is 
the  same  as  stating  that,  in  spite  of  its  obvious  shortcomings, 
the  direct  measurement  of  the  CPU  time  is  still  the  more 
reliable  way  of  comparing  different  minimization  algorithms. 

Key  Words.  Numerical  analysis,  numerical  methods,  computing 
methods,  computing  techniques,  complexity  of  computation, 
philosophy  of  computation,  comparison  of  algorithms,  computa- 
tional speed,  measurement  of  computational  speed,  number  of 
function  evaluations,  equivalent  number  of  function  evalua- 
tions, unconstrained  minimization,  mathematical  programming. 
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4.  MIELE,  A.,  BONARDO,  F.  , and  GONZALEZ,  S.,  Modifications 
and  Alternatives  to  the  Cubic  Interpolation  Process  for 
One-Dimensional  Search,  Rice  University,  Aero-Astronautics 
Report  No.  135,  1976. 

Abstract.  In  this  paper,  the  numerical  solution  of  the 
problem  of  minimizing  a unimodal  function  f (a)  is  considered, 
where  a is  scalar.  Two  modifications  of  the  cubic  interpola- 
tion process  are  presented,  so  as  to  improve  the  robustness  of 
the  method  and  force  the  process  to  converge  in  a reasonable 
number  of  iterations,  even  in  pathological  cases.  Modifica- 
tion Ml  includes  the  nonoptional  bisection  of  the  interval 
of  interpolation  at  each  iteration  of  the  process.  Modifica- 
tion M2  includes  the  optional  bisection  of  the  interval  of 

A 

interpolation:  this  depends  on  whether  the  slopes  fa(aQ)  and 

f (a  ) at  the  terminal  points  a and  a of  two  consecutive 
a o o o 

iterations  have  the  same  sign  or  opposite  sign. 

An  alternative  to  the  cubic  interpolation  process  is 
also  presented.  Th-'s  is  a Lagrange  interpolation  scheme  in 
which  the  quadratic  approximation  to  the  derivative  of  the 
function  is  considered.  The  coefficients  of  the  quadratic 
are  determined  from  the  values  of  the  slope  at  three  points: 
a^,  a2  , and  a3  = (a1  + a2)/2,  where  and  a2  are  the  end- 


points of  the  interval  of  interpolation.  The  proposed  alter- 
native is  investigated  in  two  versions,  Version  A1  and 
Version  A2.  They  differ  in  the  way  in  which  the  next  interval 
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of  interpolation  is  chosen;  for  Version  Al,  the  choice 
depends  on  the  sign  of  the  slope  fa(aQ)  f for  Version  A2,  the 
choice  depends  on  the  signs  of  the  slopes  f (a  ) and  f^ta^). 

Twenty-nine  numerical  examples  are  presented.  The 
numerical  results  show  that  both  modifications  of  the  cubic 
interpolation  process  improve  the  robustness  of  the  process. 
They  also  show  the  promising  characteristics  of  Version  A2  of 
the  proposed  alternative.  Therefore,  the  one-dimensional 
search  schemes  described  here  have  potential  interest  for 
those  minimization  algorithms  which  depend  critically  on  the 
precise  selection  of  the  stepsize,  namely,  conjugate  gradient 
methods. 

Key  Words.  One-dimensional  search,  cubic  interpolation  pro- 
cess, quadratic  interpolation  process,  Lagrange  interpolation 
scheme,  modifications  of  the  cubic  interpolation  process, 
alternatives  to  the  cubic  interpolation  process,  bisection 
process,  mathematical  programming,  interval  of  interpolation, 
numerical  analysis,  numerical  methods,  computing  methods, 
computing  techniques. 


5 . MIELE , A . , Some  Philosophical  Views  on  Algorithms  and 
Computing  Methods  in  Applied  Mathematics,  Rice  Univer- 
sity, Aero-Astronautics  Report  No.  136,  1976. 

Abstract.  This  paper  summarizes  some  of  the  work  done 
by  the  Aero-Astronautics  Group  of  Rice  University  in  the  area 
of  numerical  methods  and  computing  methods.  It  describes 
some  of  the  philosophical  thoughts  that  have  guided  this  work 
throughout  the  years.  Recommendations  are  offered  concerning 
allocation  of  funds  and  distribution  of  funds.  Additional 
recommendations  are  offered  in  order  to  bridge  the  gap  between  the 
top  management  of  government  agencies  and  the  academic  commu- 
nity. 

Key  Words.  Aerospace  engineering,  applied  mathematics, 
numerical  methods,  algorithm  research,  algorithm  development. 
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VI.  Abstracts  of  Articles  Completed 

6.  MIELE,  A.,  MOHANTY,  B.P.,  and  WU,  A.K.,  Conversion  of 
Optimal  Control  Problems  with  Free  Initial  State  into 
Optimal  Control  Problems  with  Fixed  Initial  State, 
Journal  of  the  Astronautical  Sciences,  Vol.  25,  No.  1, 
1977. 

Abstract.  This  note  considers  optimal  control  problems 
involving  the  minimization  of  a functional  subject  to  differ- 
ential constraints,  nondifferential  constraints,  initial  con- 
ditions, and  final  conditions.  The  initial  conditions  can 
be  partly  fixed  and  partly  free.  Transformation  techniques 
are  suggested,  by  means  of  which  problems  with  free  initial 
state  are  converted  into  problems  with  fixed  initial  state. 
Thereby,  it  becomes  possible  to  employ,  without  change,  some 
of  the  gradient  algorithms  already  developed  for  optimal  con- 
trol problems  with  fixed  initial  state  (for  instance,  the 
sequential  gradient-restoration  algorithm) . 

The  transformations  introduced  are  two:  (i)  a linear 
transformation  and  (ii)  a nonlinear  transformation.  In  the 
linear-quadratic  case,  the  former  preserves  unchanged  the 
basic  structure  of  the  optimization  problem,  while  this  is 
not  the  case  with  the  latter. 
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The  application  of  these  transformations  to  a problem 
of  interest  in  the  aerodynamics  of  a nonslender,  axisymmetric 
body  in  Newtonian  hypersonic  flow  is  shown.  It  consists  of 
minimizing  the  pressure  drag  for  given  values  of  the  length 
and  the  volume,  with  the  nose  radius  and  the  base  radius  being 
free.  After  transformations  (i)  and  (ii)  are  introduced,  this 
problem  is  solved  by  means  of  the  sequential  ordinary  gradient- 
restoration  algorithm  (SOGRA)  and  the  sequential  conjugate 
gradient-restoration  algorithm  (SCGRA) . 

Key  Words.  Optimal  control,  numerical  methods,  computing 
methods,  gradient-restoration  algorithms,  problems  with  free 
initial  state,  transformation  techniques,  sequential  gradient- 
restoration  algorithms,  applied  aerodynamics,  Newtonian 
hypersonic  flow,  minimum  drag  bodies. 
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7.  MIELE,  A.,  and  GONZALEZ,  S.,  On  the  Comparative 

Evaluation  of  Algorithms  for  Mathematical  Programming 
Problems,  Nonlinear  Programming  3,  Edited  by  O.L. 
Mangasarian,  R.R.  Meyer,  and  S.M.  Robinson,  Academic 
Press,  New  York,  New  York,  1978. 

Abstract.  This  paper  considers  the  comparative  eval- 
uation of  algorithms  for  mathematical  programming  problems. 

It  is  concerned  with  the  measurement  of  computational  speed 
and  examines  critically  the  concept  of  equivalent  number  of 
function  evaluations  Nq.  Does  this  quantity  constitute  a 
fair  way  of  comparing  different  algorithms? 

The  answer  to  the  above  question  depends  strongly  on 
whether  or  not  analytical  expressions  for  the  components  of 
the  gradient  and  the  elements  of  the  Hessian  matrix  are 
available.  It  also  depends  on  the  relative  importance  of 
the  computational  effort  associated  with  algorithmic  opera- 
tions vis-a-vis  the  computational  effort  associated  with 
function  evaluations. 

Both  theoretical  considerations  and  extensive  numerical 
examples  carried  out  in  conjunction  with  the  Fletcher-Reeves 
algorithm,  the  Davidon-Fletcher-Powell  algorithm,  and  the 
quasilinearization  algorithm  suggest  the  following:  the 
Ng  concept,  while  accurate  in  some  cases,  has  drawbacks  in 


other  cases?  indeed,  it  might  lead  to  a distorted  view  of 
the  relative  importance  of  an  algorithm  with  respect  to 
another . 

The  above  distortion  can  be  corrected  through  the  in- 
troduction of  a more  general  parameter,  the  time-equivalent 
number  of  function  evaluations  Ne  = T/Xg,  where  T denotes  the 
CPU  time  required  to  solve  a particular  problem  on  a parti- 
cular computer  and  Tq  denotes  the  CPU  time  required  to 
evaluate  the  objective  function  once  on  that  computer.  This 
generalized  parameter  is  constructed  so  as  to  reflect  accu- 
rately the  computational  effort  associated  with  function 
evaluations  and  algorithmic  operations. 

From  the  analyses  performed  and  the  results  obtained, 
it  is  inferred  that,  due  to  the  weaknesses  of  the  Ng  concept, 
the  use  of  the  Nq  concept  is  advisable.  In  effect,  this  is 
the  same  as  stating  that,  in  spite  of  its  obvious  shortcomings, 
the  direct  measurement  of  the  CPU  time  is  still  the  more  re- 
liable way  of  comparing  different  minimization  algorithms. 

Key  Words.  Numerical  analysis,  numerical  methods,  computing 
methods,  computing  techniques,  complexity  of  computation, 
philosophy  of  computation,  comparison  of  algorithms,  computa- 
tional speed,  measurement  of  computational  speed,  number  of 
function  evaluations,  equivalent  number  of  function  evalua- 
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8.  MIELE,  A.,  Some  Philosophical  Views  on  Algorithms  and 
Computing  Methods  in  Applied  Mathematics,  Proceedings 
of  the  Workshop  on  Decision  Information  for  Tactical 
Command  and  Control,  Airlie,  Virginia,  1976;  Edited  by 
R.M.  Thrall,  C.P.  Tsokos,  and  J.C.  Turner;  Robert  M. 
Thrall  and  Associates,  Houston,  Texas,  1976. 

Abstract.  This  paper  summarizes  some  of  the  work  done 
by  the  Aero-Astronautics  Group  of  Rice  University  in  the 
area  of  numerical  methods  and  computing  methods.  It  describes 
some  of  the  philosophical  thoughts  that  have  guided  this  work 
throughout  the  years.  Recommendations  are  offered  concerning 
allocation  of  funds  and  distribution  of  funds.  Additional 
recommendations  are  offered  in  order  to  bridge  the  gap  between 
the  top  management  of  government  agencies  and  the  academic 
community. 

Key  Words.  Aerospace  engineering,  applied  mathematics,  nume- 
rical methods,  algorithm  research,  algorithm  development. 
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9.  CLOUTIER,  J.R.,  MOHANTY,  B.P.,  and  MIELE,  A.,  Sequential 
Conjugate  Gradient-Restoration  Algorithm  for  Optimal 
Control  Problems  with  Nondifferential  Constraints,  In- 
ternational Journal  of  Control  (to  appear) . 

Abstract.  A sequential  conjugate  gradient-restoration 
algorithm  is  developed  in  order  to  solve  optimal  control 
problems  involving  a functional  subject  to  differential  con- 
straints, nondifferential  constraints,  and  terminal  constraints. 
The  algorithm  is  composed  of  a sequence  of  cycles,  each  cycle 
consisting  of  two  phases,  a conjugate  gradient  phase  and  a 
restoration  phase. 

The  conjugate  gradient  phase  involves  a single  iteration 
and  is  designed  to  decrease  the  value  of  the  functional  while 
satisfying  the  constraints  to  first  order.  During  this  iter- 
ation, the  first  variation  of  the  functional  is  minimized, 
subject  to  the  linearized  constraints.  The  minimization  is 
performed  over  the  class  of  variations  of  the  control  and  the 
parameter  which  are  equidistant  from  some  constant  multiple  of 
the  corresponding  variations  of  the  previous  conjugate  gradient 
phase.  For  the  special  case  of  a quadratic  functional  subject 
to  linear  constraints,  various  orthogonality  and  conjugacy 
conditions  hold. 

The  restoration  phase  involves  one  or  more  iterations 
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and  is  designed  to  restore  the  constraints  to  a predetermined 
accuracy,  while  the  norm  of  the  variations  of  the  control  and 
the  parameter  is  minimized,  subject  to  the  linearized  con- 
straints. The  restoration  phase  is  terminated  whenever  the 
norm  of  the  constraint  error  is  less  than  some  predetermined 
tolerance. 

The  sequential  conjugate  gradient-restoration  algorithm 
is  characterized  by  two  main  properties.  First,  at  the  end 
of  each  conjugate  gradient-restoration  cycle,  the  trajectory 
satisfies  the  constraints  to  a given  accuracy;  thus,  a se- 
quence of  feasible  suboptimal  solutions  is  produced.  Second, 
the  conjugate  gradient  stepsize  and  the  restoration  stepjsize 
can  be  chosen  so  that  the  restoration  phase  preserves  the 
descent  property  of  the  conjugate  gradient  phase;  thus,  the 
value  of  the  functional  at  the  end  of  any  cycle  is  "smaller 
than  the  value  of  the  functional  at  the  beginning  of  that  cycle. 
Of  course,  restarting  the  algorithm  might  be  occasionally  ne- 
cessary. 

To  facilitate  numerical  integrations,  the  interval  of 
integration  is  normalized  to  unit  length.  Variable- time 
terminal  conditions  are  transformed  into  fixed-time  terminal 
conditions.  Then,  the  actual  time  at  which  the  terminal  bound- 
ary is  reached  becomes  a component  of  a vector  parameter  being 
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Several  numerical  examples  are  presented,  some  pertain- 
ing to  a quadratic  functional  subject  to  linear  constraints 
and  some  pertaining  to  a nonquadratic  functional  subject  to 
nonlinear  constraints.  These  examples  demonstrate  the  feasi- 
bility as  well  as  the  rapid  convergence  characteristics  of 
the  sequential  conjugate  gradient-restoration  algorithm. 

Key  Words.  Optimal  control,  gradient  methods,  conjugate- 
gradient  methods,  numerical  methods,  computing  methods, 
gradient-restoration  algorithms,  sequential  gradient-restoration 
algorithms,  sequential  conjugate  gradient-restoration  algorithms, 
nondifferential  constraints. 


Abstract.  This  paper  considers  optimal  control  problems 
involving  the  minimization  of  a functional  subject  to  differ- 
ential constraints,  terminal  constraints,  and  a state  inequa- 


lity constraint.  The  state  inequality  constraint  is  of  a 
special  type,  namely,  it  is  linear  in  some  or  all  of  the  com- 
ponents of  the  state  vector. 

A transformation  technique  is  suggested,  by  means  of 
which  the  inequality  constrained  problem  is  converted  into  an 
equality  constrained  problem  involving  differential  constraints, 
terminal  constraints,  and  a control  equality  constraint.  The 
transformation  technique  takes  advantage  of  the  linearity  of 
the  state  inequality  constraint  so  as  to  yield  a transformed 
problem  characterized  by  a new  state  vector  of  minimal  size. 

This  concept  is  important  computationally,  in  that  the  com- 
puter time  per  iteration  increases  with  the  square  of  the 
dimension  of  the  state  vector. 

In  order  to  illustrate  the  advantages  of  the  new  trans- 
formation technique,  several  numerical  examples  are  solved  by 
means  of  the  sequential  gradient-restoration  algorithm  for 
optimal  control  problems  involving  nondifferential  constraints. 
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The  examples  show  the  substantial  savings  in  computer  time 
for  convergence,  which  are  associated  with  the  new  transforma- 
tion technique. 

Key  Words.  Optimal  control,  numerical  methods,  computing 
methods,  transformation  techniques,  sequential  gradient- 
restoration  algorithm,  nondifferential  constraints,  state 
inequality  constraints,  linear  state  inequality  constraints. 
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11.  MIELE,  A.,  MANGIAVACCHI , A.,  MOHANTY,  B.P.,  and  WU,  A.K., 
Numerical  Determination  of  Minimum  Mass  Structures  with 
Specified  Natural  Frequencies,  Rice  University,  Internal 
Memorandum,  1977. 

Abstract.  The  problem  of  the  axial  vibration  of  a can- 
tilever beam  is  investigated  both  analytically  and  numerically. 
The  mass  distribution  that  minimizes  the  total  mass  for  a given 
value  of  the  frequency  parameter  3 is  determined  using  both  the 
sequential  gradient-restoration  algorithm  (SGRA)  and  the  modi- 
fied quasilinearization  algorithm  (MQA) . Concerning  the  mini- 
mum value  of  the  mass,  SGRA  leads  to  a numerical  solution 
precise  to  5 significant  digits  and  MQA  leads  to  a numerical 
solution  precise  to  7 significant  digits. 

Comparison  of  the  optimal  beam  (a  variable-section  beam) 
with  a reference  beam  (a  constant-section  beam)  shows  that  the 
weight  reduction  associated  with  the  optimal  design  depends 
strongly  on  the  value  of  the  frequency  parameter  8.  This 
weight  reduction  is  negligible  for  B + O,  is  11.3%  for  8=1, 
is  78.6%  for  3 = 1.5,  and  approaches  100%  for  8 -*■  tt/ 2. 

Key  Words.  Structural  optimization,  dynamic  optimization, 
optimal  structures,  cantilever  beams,  sequential  gradient- 

restoration  algorithm,  modif ied-quasilinearization  algorithm, 
numerical  methods,  computing  methods. 
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12.  MIELE,  A.,  and  MOHANTY,  B.P.,  Minimax  Problems  and  the 
Computation  of  Chebyshev  Optimal  Controls , Rice  Univer- 
sity, Internal  Memorandum,  1977. 

Abstract.  A transformation  technique  has  been  developed 
in  order  to  handle  a class  of  optimal  control  problems  with  a 
Chebyshev  minimax  performance  index  within  the  framework  of 
the  classical  Mayer  problem,  wherein  the  performance  index 
becomes  an  additional  parameter  of  the  problem.  The  transfor- 
mation also  introduces  an  inequality  constraint,  which  is  akin 
to  a state  variable  inequality  constraint. 

Two  different  numerical  optimization  schemes  have  been 
explored.  The  first  one  is  an  indirect  method  involving  a 
continuous  control  and  a single  subarc.  The  second  one  is  a 
direct  method  involving  multiple  subarcs.  A two-subarc  class 
of  algorithms  has  been  developed  so  as  to  handle  minimax  prob- 
lems having  a unique  minimax  point.  Also,  a three-subarc  class 
of  algorithms  has  been  developed  so  as  to  handle  minimax  prob- 
lems having  a flat  minimax  solution:  here,  the  middle  subarc 
becomes  the  flat  minimax  solution  of  the  problem. 

The  algorithms  developed  belong  to  the  class  of  sequen- 
tial gradient-restoration  algorithms,  which  is  made  up  of  a 
sequence  of  two-phase  processes  or  cycles,  each  consisting  of 
a gradient  phase  and  a restoration  phase.  The  principal  pro- 


perty  of  the  class  of  algorithms  is  that  it  produces  a se- 
quence of  feasible  suboptimal  solutions.  To  facilitate 
numerical  implementation,  the  intervals  of  integration  have 
been  normalized  to  unit  length. 

To  illustrate  the  transformation  technique  and  the  al- 
gorithmic studies,  several  numerical  examples  are  being 
developed . 

Key  Words.  Optimal  control,  Chebyshev  optimal  control, 
numerical  methods,  computing  methods,  minimax  problems,  trans- 
formation techniques,  sequential  gradient-restoration  algorithm, 
state  inequality  constraints. 
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13.  MIELE,  A. , Transformation  Techniques  for  Optimal  Control 

Problems,  Rice  University,  Internal  Memorandum,  1977. 

Abstract.  This  paper  presents  an  analysis  of  the 
transformation  techniques  useful  in  order  to  bring  a wide 
variety  of  problems  of  optimal  control  within  the  frame  of 
three  standard  problems:  (i)  the  minimization  of  a functional 
subject  to  differential  constraints  and  terminal  constraints; 
(ii)  the  minimization  of  a functional  subject  to  differential 
constraints,  nondifferential  constraints,  and  terminal  con- 
straints; and  (iii)  the  minimization  of  a functional  subject 
to  differential  constraints,  nondifferential  constraints,  and 
multipoint  constraints. 

The  topics  treated  include:  (a)  time  normalization; 

(b)  problems  with  free  initial  state;  (c)  problems  with  bounded 
control;  (d)  problems  with  bounded  state;  and  (c)  problems  with 
bounded  time  rate  of  change  of  the  state. 

Key  Words.  Optimal  control,  numerical  methods,  computing 
methods,  gradient-restoration  algorithms,  transformation 
techniques,  time  normalization,  free  initial  state,  bounded 
control,  bounded  state,  bounded  time  rate  of  change  of  the 
state,  nondifferential  constraints,  multiple  subarcs. 


